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We analyze the constraining power of future dark energy experiments for Pseudo-Nambu Gold- 
stone Boson (PNGB) quintessence. Following the Dark Energy Task Force methodology, we forecast 
data for three experimental "stages" : Stage 2 represents in-progress projects relevant to dark energy; 
Stage 3 refers to medium sized experiments; Stage 4 comprises larger projects. We determine the 
posterior probability distribution for the parameters of the PNGB model using Markov Chain Monte 
Carlo analysis. Utilizing data generated on a ACDM cosmology, we find that the relative power 
of the different data stages on PNGB quintessence is roughly comparable to the DETF results for 
the wo — w a parametrization of dark energy. We also generate data based on a PNGB cosmological 
model that is consistent with a ACDM fiducial model at Stage 2. We find that Stage 4 data based 
on this PNGB fiducial model will rule out a cosmological constant by at least 3a. 



I. INTRODUCTION 

A growing number of observations indicate that the ex- 
pansion of the universe is accelerating. Given our current 
understanding of physics, this phenomenon is a mystery. 
If Einstein's gravity is correct, then it appears that ap- 
proximately 70 percent of the energy density of the uni- 
verse is in the form of a "dark energy" . Although there 
are many ideas of what the dark energy could be, as of 
yet, none of them stands out as being particularly com- 
pelling. Future observations will be crucial to developing 
a theoretical understanding of dark energy. 

A number of new observational efforts have been pro- 
posed to probe the nature of dark energy, but evaluating 
the impact of a given proposal is complicated by our poor 
theoretical understanding of dark energy. In light of this 
issue, a number of model independent methods have been 
used to explore these issues (see for example [H, 0). No- 
tably, the Dark Energy Task Force (DETF) produced a 
report in which they used the w — w a parametrization 
of the equation of state evolution in terms of the scale 
factor, w(a) = w + u? a (l — a) ||. The constraints on 
the parameters w and w a were interpreted in terms of a 
"figure of merit" (FoM) designed to quantify the power 
of future data and guide the selection and planning of 
different observational programs. Improvements in the 
DETF FoM between experiments indicate increased sen- 
sitivity to possible dynamical evolution of the dark en- 
ergy. This is crucial information, since current data is 
consistent with a cosmological constant and any detec- 
tion of a deviation from a cosmological constant would 
have a tremendous impact. There are, however, a num- 
ber of questions left unanswered when the dark energy is 
modeled using abstract parameters such as wg — w a that 
are perhaps better addressed with the analysis of actual 
theoretical models of dark energy. 

First of all, the w — w a parametrization is simplistic 
and not based on a physically motivated model of dark 
energy. Although simplicity is part of the appeal of this 
parametrization, some of the most popular dark energy 
models exhibit behavior that cannot be described by the 



w —w a parametrization. The PNGB Quintessence model 
considered in this paper, for instance, allows equations 
of state that cannot be approximated by the w — w a 
parametrization. 
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FIG. 1: Examples of possible equations of state for PNGB 
Quintessence (left panel). Attempts to imitate this behavior 
with w(z) curves for the wo — w a parametrization are depicted 
on the right. 

Conversely, the w — w a parametrization may allow so- 
lutions that do not correspond to a physically motivated 
model of dark energy. Because of these issues, one could 
wonder whether the DETF FoM's are somehow mislead- 
ing. Various concerns with the DETF analysis have al- 
ready been explored. Albrecht and Bernstein used a more 
complex parametrization of the scale factor to check the 
validity of the u>o — w a approximation and Huterer 
and Peiris considered a generalized class of scalar field 
models [5J . Additionally, it has been suggested that the 
DETF data models might be improved (see for instance 
0). 

As of yet, however, no actual proposed models of dy- 
namical dark energy have been considered in terms of 
future data. Given the issues above, such an analysis is 
an important compliment to existing work. Of course 
all specific models of dark energy are suspect for various 
reasons, and one can just as well argue that it is better 
to make the case for new experiments in a more abstract 
parameter space rather than tying our future efforts to 
specific models that themselves are problematic. Rather 
than "take a side" in this discussion, our position is that 
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given the diversity of views on the subject, a model-based 
assessment will have an important role in an overall as- 
sessment of an observational program. 

In this paper we consider the pseudo-Nambu Goldstone 
boson quintessence model of dark energy Q. As one of 
the most well-motivated quintessence models from a par- 
ticle physics perspective, it is a worthwhile one to study. 
We use the data models forecasted by the DETF and gen- 
erate two types of data sets, one based on a AC DM back- 
ground cosmology and one based on a background cos- 
mology with PNGB dark energy using a specific fiducial 
set of PNGB parameters. We determine the probabil- 
ity space for the PNGB parameters given the data using 
Markov Chain Monte Carlo analysis. This paper is part 
of a series of papers in which a number of quintessence 
models are analyzed in this manner [8, 9]. 

We show that the allowed regions of parameter space 
shrink as we progress from Stage 2 to Stage 3 to Stage 4 
data in much the same manner as was seen by the DETF 
in the u>o — w a space. This result holds for both KCDM 
and the PNGB data models. Additionally, with our 
choice of PNGB fiducial background model, we demon- 
strate the ability of Stage 4 data to discriminate between 
a universe described by a cosmological constant and one 
containing an evolving PNGB field. As cosmological data 
continues to improve, careful analysis of specific dark en- 
ergy models using real data will become more and more 
relevant. MCMC analysis can be computationally inten- 
sive and time-consuming. Since future work in this area 
is likely to encounter similar challenges, we discuss some 
of the difficulties we discovered and solutions we imple- 
mented in our MCMC exploration of PNGB parameter 
space. 



II. PNGB QUINTESSENCE 

Quintessence models of dark energy are popular con- 
tenders for explaining the current acceleration of the uni- 
verse 

nam 

12]. Although the cosmological constant is 
regarded by many to be the simplest theory of the dark 
energy, the required value of the cosmological constant 
appears to be many orders of magnitude too small in 
naive particle theory estimates. In quintessence mod- 
els this problem is not solved. Instead it is generally 
sidestepped by assuming some unknown mechanism sets 
the vacuum energy to exactly zero, and the dark energy 
is due to a scalar field evolving in a pressure dominated 
state. As such fields can appear in many proposed "fun- 
damental theories", and as the mechanism mimics ideas 
familiar from cosmic inflation, quintessence models are 
regarded by many (but certainly not by everyone fl3j |) 
to be at least as plausible as a cosmological constant 
Pill- 

Here the quintessence field is presumed to be homoge- 
neous in space, and is described by some scalar degree of 
freedom <f> and a potential V(<p) which governs the field's 
evolution. In an FRW spacetime, the field's evolution is 



given by 

^ + 3^0+^ = (1) 

where 



and 

H 2 = J^jliPr + Pm + P<f, + Pk) (3) 

where Mp is the reduced Planck mass, p r is the energy 
density of radiation, p m is the energy density of non- 
relativistic matter and pk is the effective energy density 
of spacetime curvature. The energy density and pressure 
associated with the field are 

P4> = \fi + V{4>), P^\tf-V{<\>) (4) 

and the equation of state w is given by 




If the potential energy dominates the energy of the field, 
then as can be seen in Eq. [?]the pressure will be negative 
and in some cases can be sufficiently so to give rise to 
acceleration as the universe expands. 

The PNGB model of quintessence is considered com- 
pelling because it is one of the few models that seems 
natural from the perspective of 4-D effective field the- 
ory. In order to fit current observations, the quintessence 
field must behave at late times approximately as a cos- 
mological constant. It must be rolling on its potential 
without too much contribution from kinetic energy, and 
the value of the potential must be close to the observed 
value of the dark energy density, which is on the order 
of the critical density of the universe p c = SH^m 2 = 
1.88 x W~ 26 h 2 kg/m 3 or 2.3 x lCT 120 /) 2 in reduce Planck 
units, where h = Ho/100. These considerations re- 
quire the field to be nearly massless and the potential 
to be extraordinarily flat from the point of view of par- 
ticle physics. In general, radiative corrections generate 
large mass renormalizations at each order of perturba- 
tion theory unless symmetries exist to suppress this ef- 
fect [HI, [I?], [l8| . In order for such fields to seem reason- 
able, at least on a technical level, their small masses must 
be protected by symmetries such that when the masses 
are set to zero they cannot be generated in any order of 
perturbation theory. Many believe that Pseudo-Nambu- 
Goldstone bosons are the simplest way to have ultra-low 
mass, spin-0 particles that are natural in a quantum field 
theory sense. An additional attraction of the model is 
that the parameters of the PNGB model might be re- 
lated to the fundamental Planck and Electroweak scales 
in a way that solves the cosmic coincidence problem [l9| . 
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The potential of the PNGB field is well-approximated 



V = Af 4 [cos(p + l] 



(6) 



(where higher derivative terms and instanton corrections 
are ignored). The evolution of the dark energy is con- 
trolled by the two parameters of the PNGB potential, 
M 4 and /, and the initial conditions, <f>i and (f>i. Wc 
take </>/ = 0, since we expect the high expansion rate of 
the early universe to rapidly damp nonzero values of 4>i ■ 
The initial value of the field, <j>i, takes values between 
and 7r/. This is because the potential is periodic and 
symmetric. Since a starting point of 4>i / f on the po- 
tential is equivalent to starting at mr — cf>j / f and rolling 
in the opposite direction down the potential, we require 

< 0/// < 7T. 
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FIG. 2: An example of a PNGB model and its resulting cos- 
mological solution for typical values of the PNGB parameters. 
The PNGB potential V((j>) (dashed curve, right panel) is in 
units of h 2 . The evolution of the PNGB field along the po- 
tential since the radiation era is shown by the solid curve 
overlaying the potential. The energy densities in the right 
panel are related to these units via Q, = uj/h 2 . 

Additionally, we place the bound / < M p . As will 
be discussed in the following section, this is necessary to 
cut off a divergent direction so that the MCMC chains 
converge. There are theoretical reasons for this bound 
as well. For one, it is valid to neglect higher derivative 
terms of the PNGB potential, Eq. ©, at least as long 
as / < M p . In general, we don't expect to understand 
4-D effective field theory at energies much larger than 
this. Additionally, there are indications from string the- 
ory that / cannot be larger than M p (20L l2lj ] . 



III. 



ANALYSIS AND MCMC 



Following the DETF methodology, we generate data 
sets for future supernova, weak gravitational lensing, 
baryon acoustic oscillation and cosmic microwave back- 
ground observations. These observations are forecasted 
for three experimental "stages": Stage 2 represents in- 
progress projects relevant to dark energy; Stage 3 refers 
to medium-sized experiments; Stage 4 comprises larger 



dark energy projects including a large ground-based sur- 
vey (LST), and/or a space-based program. (Stage I, not 
considered in our analysis, represents already completed 
dark energy experiments, and is less constraining than 
Stage 2 data.) "Optimistic" and "pessimistic" versions 
of the same simulated data sets give different estimates 
of the systematic errors. More information on the spe- 
cific data models is given in Appendix [5] (or see also the 
technical appendix of the DETF report Q). In our work 
we did not use the cluster data because of the difficulty of 
adapting the DETF construction to a quintessence cos- 
mology (the same reasons given in Q). 

We generate and analyze two versions of data. One is 
built around a cosmological constant model of the uni- 
verse. The other is based on a PNGB fiducial model. 
The latter is chosen to be consistent with a cosmological 
constant for Stage 2 data. 

We use Markov Chain Monte Carlo a naly sis with a 
Metropolis-Hastings stepping algorithm j22l [23l . [24| to 
evaluate the likelihood function for the parameters of our 
model. The details of our methods are discussed in Ap- 
pendix |Bj MCMC lends itself to our analysis because 
our probability space is both non-Gaussian and also de- 
pends on a large number of parameters. These include 
the PNGB model parameters: M 4 , /, and <j>i : the cos- 
mological parameters: w m , ujk, wb, 5^, n s (as defined 
by the DETF), and the various nuisance and/or photo-z 
parameters accounting for error and uncertainties in the 
data). 

In order for the results of an MCMC chain to be mean- 
ingful, there must exist a finite, stationary distribution to 
which the Markov chain may converge in a finite number 
of steps. Degeneracies between parameters, i.e. combi- 
nations of different parameters that give rise to identical 
cosmologies, correspond to unconstrained directions in 
the probability distribution. Unless some transformation 
of parameters is found and/or a cut-off placed on these 
parameters, the MCMC will step infinitely in this direc- 
tion and can never converge to a stationary distribution. 
Additionally, the shape of the probability distribution 
can drastically effect the efficiency of the chain. A large 
portion of the task of analyzing the PNGB model, there- 
fore, involves finding convenient parameterizations and 
cutoffs to facilitate MCMC exploration of the posterior 
distribution. 

The probability space of the PNGB model becomes 
more tractable from an MCMC standpoint if we trans- 
form from the original variables to ones that are more 
directly related to cosmological observables constrained 
by the data. Such parameterizations make it easier to 
identify degeneracies and also tend to make the shape of 
the probability distribution more Gaussian. As discussed 
in Section [Til the dynamics of the PNGB field depend on 
its potential V{4>) = M i (cos(j) + f), and the specific 

values of M , / and 0/. In order to fit current data the 
field must hang on the potential approximating a cosmo- 
logical constant for most of the expansion history of the 
universe. If the field never rolls it acts as a cosmological 
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constant for all times with a value corresponding to the 
initial energy density of the field, Vj — V(cj)j). To first 
order, then, Vj sets the overall scale of the dark energy 
density. Since Vi has more physical significance than M 4 , 
it is a more efficient choice for our MCMC analysis. 

Additionally, the "phase", 4>i/f, of the field's starting 
point in the cosine potential is closely related to the ini- 
tial slope of the potential. The slope affects the timescale 
on which the field will evolve. If, for instance, <fii / f = 
the field starts out on the very top of the potential where 
the slope is exactly zero, and the field will not evolve. 
Starting closer to the inflection point of the potential 
results in a steeper initial slope and the field will roll 
faster toward the minimum. Since the variable 4>i can 
correspond to both flat initial slopes (if / is large) or 
steep ones (if / is small), 4>i / f is more directly related 
to the dynamics of the field and is therefore a superior 
parameter choice. These new parameters, as illustrated 
in Fig. [3j also result in a probability distribution that is 
more Gaussian and thus is more easily explored by our 
MCMC algorithm. 
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FIG. 4: PNGB potentials (dashed) and with the entire field 
evolution shown in thick solid curves. The different curves 
show increasing values of / from left to right. The smallest 
value of / (bottom curve) gives a nearly static dark energy 
and fits a cosmological constant well. Any larger value for / 
will also fit the data because the potential will be flatter, and 
the field will evolve even less. 
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FIG. 3: 2-D confidence regions for / vs. M 4 (left panel) 
and / vs. Vi (right panel). The concave feature of the / 
vs. M 4 contours means it is inefficiently explored by MCMC. 
Contours in the / — Vi space are nearly Gaussian and better 
facilitate convergence. 

Even more important than choosing physically rele- 
vant parameters, is deciding how to handle divergent di- 
rections in probability space. For PNGB quintessence 
the parameter / must be cut off in some way because it 
can become arbitrarily large without being constrained 
by the data. For the fiducial model based on a ACDM 
universe, solutions where the field does not evolve for the 
entire expansion history of the universe, i.e. behaves as a 
cosmological constant, can fit the forecast data perfectly. 
If such a choice of parameters is found, then larger and 
larger values of / will only make the potential flatter and 
flatter. If the field did not evolve significantly for the 
smaller values of /, this will be even more true as the 
potential flattens. Hence, / can become arbitrarily large 
and the cosmological observables will remain identical. 

In general, it is possible to achieve identical deviations 
from a cosmological constant by increasing / while at 
the same time moving </>/// towards the inflection point 
of the potential. Increasing / flattens the potential, but 
by changing <fii/f, the slope of the potential can be held 



nearly constant and the evolution of the field will not 
change significantly. In order to achieve results with 
MCMC, it is necessary to choose some cutoff / so that 
this infinite direction is bounded. We choose / < M p be- 
cause there is some theoretical motivation for this choice 
as detailed in Sec. [TTJ 



IV. RESULTS 

A. ACDM Fiducial Model 

In this section we present the results of our MCMC 
analysis for the combined data sets based on a ACDM 
fiducial model. The model parameters (represented in 
Table []} are in units of h 2 . 



TABLE I: ACDM Fiducial Parameter Values (energy densi- 
ties in units of h 2 ) 
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Stage 2 combines supernovae, weak lensing and CMB 
data. Stages 3 and 4 additionally include BAO data. We 
marginalize over all but two parameters to calculate 2-D 
contours for parameters of the PNGB model and find the 
68.27%, 95.44% and 99.73% (1, 2 and 3 sigma) confidence 
regions. 

Fig. [5] depicts the contours in the Vj — cf>ij f plane for 
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Stage 2, and the optimistic versions of Stage 3, Stage 4 
space and Stage 4 LST-ground combined data. The hori- 
zontal axis where f = corresponds to a cosmological 
constant. (As explained above, the field is starting ex- 
actly at the top of its potential and does not roll because 
the potential is flat at this point.) 
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FIG. 5: Vi — 4>i/f 1, 2 and 3 sigma confidence regions for 
DETF "optimistic" combined data sets. 
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FIG. 6: / — 4>i/f 1, 2 and 3 sigma confidence regions for 
DETF "optimistic" combined data sets. 

The value of Vj on this axis, therefore, represents the 
dark energy density, ujde, or A. The contours, as ex- 
pected, are centered around Vj = .38, the fiducial value 
of lode- ft can be seen that the area of the contours 
shrinks from Stage 2 to Stage 3 and again from Stage 3 
to Stage 4. The shrinking in the <fii/f direction roughly 
corresponds to constraining deviations from a cosmolog- 
ical constant. (Although this interpretation is a slight 
oversimplification, since for larger values of /, (j)i/f can 
be non-zero and perceptible deviations from a cosmologi- 
cal constant will not occur until sometime in the future.) 
The reduction in the Vj direction reflects constraints the 



data places on the contribution from the dark energy to 
the energy density of the universe. 

Figs. [6] depicts the f — <t>i/f contours. Although all 
values of / are allowed, as / approaches zero the PNGB 
potential gets narrower, and the phase must start closer 
to zero, or else the field will evolve too quickly to its 
vacuum state. (The very thin part of the distribution 
close to / = is not resolved by the MCMC analysis.) 
For larger values of /, 4>i/f may start further from the 
peak of the potential without the field evolving much. 
Even for Stage 2, however, the field may not start past 
the inflection point of the potential. 

Often it is assumed that the PNGB field is initially 
displaced a small amount from the potential minimum. 
But with the constraint we have placed on /, this region 
of parameter space is no longer accessible for data based 
on a KCDM fiducial model. This is because as the field 
starts lower down the slope of the potential, the peak of 
the potential must be raised so that Vj may reflect the ap- 
proximate energy density needed by the dark energy. But 
as the peak of the potential gets higher it also becomes 
steeper (since / is bounded) and the field evolves too 
quickly to fit the data. It has been suggested, however, 
that since we expect quantum fluctuations to displace 
the field from the top of potential, it is more reasonable 
to expect the PNGB field to start after the inflection 
point of the potential 16]. If either this argument or 
the theoretical reasons for the bound / < M p could be 
made more convincing, experimental results consistent 
with a cosmological constant could potentially rule out 
the PNGB model. As it stands, however, we do not feel 
the arguments constraining / and <pi are robust enough 
to make such a claim. 
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FIG. 7: Vi — Sujde 1, 2 and 3 sigma confidence regions for 
optimistic combined data. Here Sujde is the amount of change 
in the dark energy density from the radiation era until today. 



Fig. [7] depicts Vj versus Slvde, where 5u>be = Vz — 
ude{o- — !)■ Since PNGB quintessence is a "thawing" 
model of dark energy [25|, that is, it starts as a cosmo- 
logical constant until the field begins to roll causing the 
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amount of dark energy to decrease, Slode reflects the 
amount the dark energy has deviated from a cosmologi- 
cal constant. As the DETF found, subsequent stages of 
data do better at constraining the evolution of the dark 
energy. The fact that Stage 4 space seems a little more 
constraining than ground reflects the fact that ground 
and space data are sensitive to slightly different features 
in the dark energy evolution and will be more or less 
powerful at different redshifts. Other quintessence mod- 
els, such as the Albrecht-Skordis model @ are somewhat 
better constrained by the DETF Stage 4 ground data 
models than by DETF Stage 4 space. 



B. PNGB Fiducial Model 

In addition to considering a AC DM fiducial model, we 
evaluate the power of future experiments assuming the 
dark energy is really due to PNGB quintessence. Our 
PNGB fiducial parameter values (shown in Table [IT] in 
units of h 2 ) were chosen such that the fiducial model lies 
within the 95% confidence region for Stage 2 KCDM 
data, but demonstrates a small amount of dark energy 
evolution that can be resolved by Stage 4 experiments. 



TABLE II: Fiducial Parameter Values (energy densities and 
Vr in units of h 2 ,f in reduced Planck units 
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FIG. 9: Vi — </>/l,2 and 3 sigma confidence regions for DETF 
optimistic combined data sets using the PNGB background 
cosmological model. 
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FIG. 10: / — (j>i/f 1, 2 and 3 sigma confidence regions for 
DETF optimistic combined data sets using the PNGB back- 
ground cosmological model. 




FIG. 8: The evolution of the PNGB fiducial model field (left 
panel, solid curve) in the PNGB potential (dashed curve). 
The corresponding equation of state evolution is shown in the 
right panel. 

The left panel of Fig. [8] shows the potential and evolu- 
tion of the field for this model. The right panel depicts 
w(z). It can be seen that today (z — 0) the deviation of 
the field from w(z) = — 1 is only about 10%. 



Repeating our MCMC analysis for the PNGB fiducial 
model, we again marginalize over all but two parameters 
to depict the 2-d confidence regions for the dark energy 
parameters. Fig. [5] and depict the Vj — <j>j/f contours. 
It can be seen that the 4>i / f = axis corresponding 
to the field sitting on the top of its potential and not 
evolving, is allowed at Stage 2 but becomes less favored 
by subsequent stages of the data. By Stage 4 it is ruled 
out by more than 3er. 

Fig. [TOldepicts the (fii/f — f contours. Again it can be 
seen that for larger values of /, </>/// must be non-zero. 
By Stage 4 optimistic, only extreme fine tuning with / 
allows 4>i I f to approach zero, so that the field will be 
displaced from the top of the potential and have started 
to roll by just the right amount by late times. 

Fig. [TT1 depict Vi versus Slude- At Stage 2 Sujde = 
is still is within the la confidence region. But subsequent 
stages of the data disfavor this result. Stage 4 optimistic 
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FIG. 11: Vj — Sujde 1, 2 and 3 sigma confidence regions for 
DETF optimistic combined data, where Sujde is the amount 
of change in the dark energy density from the radiation era 
until today. 



rules out zero evolution of the dark energy by more than 
3cr. 



V. DISCUSSION AND CONCLUSIONS 

With experiments such as the ones considered by the 
DETF on the horizon, data sets will be precise enough 
to make it both feasible and important to analyze dy- 
namic models of dark energy. The analysis of such mod- 
els, therefore, should play a role in the planning of these 
experiments. 

With our analysis of PNGB quintessence, we have 
shown how future data can constrain the parameter space 
of this model. We have shown likelihood contours for a 
selection of combined DETF data models, and found the 
increase in parameter constraints with increasing data 
quality to be broadly consistent with the DETF results 
in wo — w a space. Direct comparison with the DETF fig- 
ures of merit is non-trivial because PNGB quintessence 
depends on three parameters, whereas the DETF FoM 
were calculated on the bases of two, but in our two di- 
mensional projections we saw changes in the area that are 
consistent with DETF results. Specifically, the DETF 
demonstrated a factor of roughly three decrease in al- 
lowed parameter area when moving from Stage 2 to good 
combinations of Stage 3 data, and a factor of about ten 
in area reduction when going to from Stage 2 to Stage 
4. We saw decreases by similar factors in our two dimen- 
sional projections. We have presented likelihood contour 
plots for specific projected data sets as an illustration. 
In the course of this work we produced many more such 
contour plots to explore the other possible data combi- 
nations considered by the DETF including the data with 
"pessimistic" estimates of systematic errors. We found 
no significant conflict between our results in the PNGB 
parameter space and those of the DETF in WQ — w a space. 



As discussed in [15j], we believe the fact that we have 
demonstrated (here and elsewhere [1,0) results that are 
broadly similar to those of the DETF despite the very 
different families of functions w(a) considered is related 
to the fact pointed out in [J] that overall the good DETF 
data sets will be able to constrain many more features of 
w(a) than are present in the wq — w a ansatz alone. 

As data continues to improve, MCMC analysis of dy- 
namic dark energy models will likely become more pop- 
ular. Our experience with the PNGB model could be 
relevant to future work. We find that the theoretical pa- 
rameters of the model are not in general the best choice 
for MCMC. Transforming to variables that are closely 
related to the physical observables can help MCMC con- 
verge more efficiently. Additionally, it is necessary to cut 
off unconstrained directions in parameter space. It would 
be desirable to find bounds that have some physical mo- 
tivation. For PNGB quintessence, we find that the initial 
value of the potential, V (</>/), and the initial "phase" of 
the field, </>///, are more convenient than the original 
model parameters, and that there is some motivation for 
placing the bound / < M v . 

Finally, we have demonstrated the power Stage 4 data 
will have for detecting time evolution of the dark energy. 
The PNGB fiducial model we choose is consistent with 
Stage 2 data (and with current data by extension). If, 
however, the universe were to in fact be described by such 
a dark energy model, then by Stage 4 we would know to 
better than 3 sigma that there is a dynamic component 
to the dark energy. 



APPENDIX A: DATA 

For each step in the MCMC chain, we integrate numer- 
ically to calculate the theoretical quantities dependent on 
the dark energy. We start our integration at early times 
with a = 10 -15 and we end the calculation at a = 2. 
We compare these value with the observables generated 
based on our fiducial models. With the uncertainties in 
the data forecast by the DETF we can calculate the like- 
lihood for each step in the chain. What follows is an 
overview of the likelihood calculation for each type of 
observation we consider. 



1. Type la Supernovae 

After light curve corrections, supernovae observations 
provide the apparent magnitudes, rrij, and the redshift 
values, Zi, for supernova events. The apparent magni- 
tudes are related to the theoretical model through the 
distance modulus, /Lt(-Zj), by 



where 



m l = M + (j,(zi) 



H(zi) = 5log w (di(zi)) + 25 



(Al) 



(A2) 
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M is the absolute magnitude and 



2. Baryon Acoustic Oscillations 



f i amh(y/\k\x{zi)) k<0 



dh(zi) = - < 



X {zi) _ k = 

L ^sin(v/j%(*i)) fc>0 



(A3) 



and 



X(^) = % - *7(z») 



da 



. a?H{a) 



(A4) 



with \k\ = H*\n k \ = (^) 2 K|. 

Uncertainties in absolute magnitude M as well as the 
absolute scale of the distance modulus lead to the intro- 
duction of an offset fj,off nuisance parameter in all SNc 
data sets, giving n{zi) — > /i(zi) + Mo//- 

Other systematic errors are modeled by more nuisance 
parameters. The peak brightness of supernovae, for in- 
stance, may have some z-dependent behavior that is not 
fully understood. We include this uncertainty in our 
analysis by allowing the addition of small linear and 
quadratic offsets in z. Additionally, each SNe data model 
combines a collection of nearby supernovae with a collec- 
tion of more distant ones. Possible differences between 
the two groups are modeled by considering the addition 
of a constant offset to the near group. The distance mod- 
ulus becomes 

f^{ z i)calc M( z i) "f" Mo// ~t~ Hlin z i ~t~ Hquad z i ~\~ Hshift z near 

(A5) 

In addition, some experiments will measure super- 
novae redshifts photometrically instead of spectroscopi- 
cally. There may be a bias in the measurement of the 
photo-z's in each bin. This uncertainty is expressed 
by another set of nuisance parameters, Szi, that can 
shift the values of each Zj. These observables become 

Mj = fl(Zi + 5Zi). 

Priors are assigned to each of the nuisance parameters 
(except for Mo//, which is left unconstrained) which re- 
flect the projected strength of the various observational 
programs. Additionally, statistical errors are presumed 
to be Gaussian and are given by the diagonal covariancc 
matrix Cij = of , where o» reflects the uncertainty in the 
Hi observables for each data set. 

The likelihood function L for the supernovae data can 
be calculated from the chi-squared statistic, where x 2 = 
—2lnL. For data sets with photometrically determined 
redshifts chi-squared is 



X 



H(Zi) - H( z t)data^ i^lin f-Lquad 

) + T2 r- — 2 



®lin quad 



5i 



linear \~"V * \ 

^2 + 2-J*—>- 



(A6) 



For data sets with spectroscopic redshifts the chi-squared 
is the same minus the contribution from redshift shift 
parameters. 



Large scale variations in the baryon density in the uni- 
verse has a signature in the matter power spectrum that 
when calibrated via the CMB provides a standard ruler 
for probing the expansion history of the universe. The 
observables for BAO data (after extraction from the mass 
power spectrum) are the comoving angular diameter dis- 
tance, d c °{zi), and the expansion rate, H(zi), where 



dl 



adi 



(A7) 



and Zi indicates the z bin for each data point. The qual- 
ity of the data probe is modeled by the covariance matrix 
for each observable type, as described in section 4 of the 
DETF technical appendix. Additionally, some BAO ob- 
servations use photometrically determined redshifts, in 
which case Szi are added as nuisance parameters as for 
the supernovae, to describe the uncertainty in each red- 
shift bin. 

The likelihood function for BAO observations is 



X 



^ f dl°(z i )-dl°_ data (z i ) H(Zj) - H da ta(Zi) ^ 
2-^ K a* + rr2 ' 



+ 



Szi 



(A8) 



3. Weak Gravitational Lensing 

Light from background sources is deflected from a 
straight path to the observer by mass in the foreground. 
From high resolution imaging of large numbers of galax- 
ies, it is possible to detect statistical correlations in the 
stretching of galaxies, "cosmic shear". From this fore- 
ground mass distributions can be determined. The mass 
distribution as a function of redshift provides a probe of 
the growth history of density perturbations, g(z), where 
g(z) (in linear perturbation theory) depends on dark en- 
ergy via 



g + 2Hg 



3ri m i? 2 

2a 3 



-9- 



(A9) 



Additionally, because the amount of lensing depends on 
the ratios of distances between the observer, the lens and 
the source, gravitational lensing also probes the expan- 
sion history of the universe, D(z). 

The direct observables of weak lensing surveys consid- 
ered by the DETF are the power spectrum of the lensing 
signal and the cross-correlation of the lensing signal with 
foreground structure. Systematic and statistical uncer- 
tainties are described by a Fisher matrix in this space. As 
is detailed in the DETF appendix, it is possible to trans- 
form from this parameter space to the variables directly 
dependent on dark energy, g(z) and D(z). These become 
the weak lensing observables we use in our analysis. 

In addition to depending on the dark energy model, 
weak lensing observations depend on the cosmological 
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parameters io m , matter density, u>b, Baryon density, Wk, 
effective curvature density, ng, the spectrai index, and 
lnA s , the amplitude of the primordial power spectrum. 
These parameters are treated as nuisance parameters 
with priors imposed by the Fisher matrix. 

Lastly, since ground-based lensing surveys will pho- 
tometrically determine redshifts, as for SNe and BAO 
data, we must model the uncertainty in redshift bins. 
Again this is done by allowing each bin to vary by 
some amount Szi. 

The weak lensing observables are given be the vector 

X obs = [(u m ,uj k ,uiB,n s ,S( : ),d L a (z i ),g(z i ),ln(a(z l )} 

(A10) 

where a(zi) is the scale factor corresponding to the red- 
shift bins for the data. The error matrix is non-diagonal 
in the space of these observables so chi-squared is given 
by 

\ obs ~~ obs — data )Fl ens i n g (~jt G bs " obs — data) 

(All) 



4. Planck CMB 

As with baryon oscillations, observations of 
anisotropics in the cosmic microwave background 
probe the expansion history of the universe by providing 
a characteristic length scale at the time of last scattering. 
As with weak lensing, our Planck observables are ex- 
trapolated from the CMB temperature and polarization 
maps. The observable space constrained by Planck 
becomes: n s , ui m , u>b, 5^, ln(&s). These variables are 
constrained via the Fisher matrix in this space (we use 
the same one used in Q). The chi-squared is calculated 
as 

(Xobs - X obs-dataV 
(A12) 



APPENDIX B: MCMC 

Markov Chain Monte Carlo simulates the likelihood 
surface for a set of parameters by sampling from the pos- 
terior distribution via a series of random draws. The 
chain steps semi-stochastically in parameter space via the 
Metropolis-Hastings algorithm such that more probable 
values of the space are stepped to more often. When the 
chain has converged it is considered a "fair sample" of 
the posterior distribution, and the density of points rep- 
resents the true likelihood surface. (Explanations of this 
technique can be found in H H H^] ) • 

With the Metropolis-Hastings algorithm, the chain 
starts at an arbitrary position 9 in parameter space. A 
candidate position 9' for the next step in the chain is 
drawn from a proposal probability density q(9,9'). The 



candidate point in parameter space is accepted and be- 
comes the next step in the chain with the probability 



a(6,9') = min{l 



P(0')q(O',O) 
P(9)q{9,9') 



} 



(Bl) 



where P{9) is the likelihood of the parameters given the 
data. If the proposal step 9' is rejected, the point 9 be- 
comes the next step in the chain. Although many dis- 
tributions are viable for the proposal density q(9,9'), for 
simplicity we have chosen to use a Gaussian normal dis- 
tribution. (It should be noted that, in general, the dark 
energy parameters of the model are not Gaussian dis- 
tributed. The power of the MCMC procedure lies in the 
fact that it can probe posterior distributions that are 
quite different from the proposal density q(9,9').) Since 
this is symmetric, q(9,9') = q(9',9), we need only con- 
sider the ratios of the posteriors in the above stepping 
criterion. 

For the results of the Markov chain to be valid, it must 
equilibrate, i.e. converge to the stationary distribution. 
If such a distribution exists, the Metropolis-Hastings al- 
gorithm guarantees that the chain will converge as the 
chain length goes to infinity. In practice, however, we 
must work with chains of finite length. Moreover, from 
the standpoint of computational efficiency, the shorter 
our chains can be and still reflect the true posterior dis- 
tribution of the parameters, the better. Hence a key 
concern is assuring that our chains have equilibrated. 
Though there are many convergence diagnostics, chains 
may only fail such tests in the case of non-equilibrium; 
none guarantee that the chain has converged [28(. We 
therefore monitor the chains in a variety of ways to con- 
vince ourselves that they actually reflect the underlying 
probability space. 

Our first check involves updating our proposal distri- 
bution q(6, 9'), which we have already chosen to be Gaus- 
sian normal. Each proposal step is drawn randomly from 
this distribution. The size of the changes generated in 
any given parameter direction depend on the covariancc 
matrix we use to define q(9,9'). We start by guessing 
the form of the covariance matrix and run a short chain 
(O(10 5 ) steps) after which we calculate the covariance 
matrix of the Markov chain. We then use this covariance 
matrix to define the Gaussian proposal distribution for 
the next chain. We repeat this process until the covari- 
ance matrix stops changing systematically. This implies 
that the Gaussian approximation to the posterior has 
been found. In addition to indicating convergence, this 
also assists the efficiency of our chains. The more the 
proposal distribution reflects the posterior, the quicker 
the Markov chain will approximate the underlying distri- 
bution. 

One convergence concern is that we might not be ex- 
ploring the entire probability space. It is possible, for 
instance, that if we started our chains at a random point 
in probability space and our step sizes are too small, the 
chain may have wandered to a local maximum from which 
it will not exit in a finite time. We could be missing other 
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features of the underlying probability space. We convince 
ourselves that this is not the case by starting chains at 
different points in parameter space. We find that the 
chains consistently reflect the same probability distribu- 
tion, and hence we conclude that we are truly sampling 
from the full posterior. 

After we have determined our chains are fully explor- 
ing probability space and we have optimized our Gaus- 
sian proposal distribution, we run a longer chain to bet- 
ter represent the probability space of our variables. We 
consider the chain to be long enough when the 95% con- 
tour is reasonably smooth. For most data sets, chains 
of O(10 6 ) are sufficient although the larger the probabil- 
ity space, the longer the chains must be. (In particular, 
Stage 4 ground data involves a large number of nuisance 
parameters and may take two to three times longer to 
return smooth contours.) With the final chains, we must 
control for both burn-in and correlations between param- 
eter steps. Burn-in refers to the number of steps a chain 
must take before it starts sampling from the stationary 
distribution. Because we have already run a number of 
preliminary chains, we know approximately the mean pa- 
rameters of our model. We find that the means refer to 
a point in probability space close to the maximum of the 
distribution. (Generically, this does not have to be true 
if the probability space is asymmetrical.) If we use this 
as our starting point, our chains do not have to wander 
long before they appear to sample from the stationary 
distribution. We control for this by removing different 
amounts from the start of the chain. For instance, if 
we cut out the first 1000 steps and calculate the con- 
tours and compare this to contours calculated with the 
first 100000 steps removed we find that the shape of the 
2 — D contours remain essentially the same. We can con- 
clude, therefore, that chains very quickly begin sampling 
the posterior distribution and we need not worry about 
burn-in. 

Correlations between steps may also effect the repre- 
sentativeness of the samples generated via MCMC. The 
effects, however, may be controlled for by either thinning 
the chains by a given amount or by running chains of 



sufficient length such that the correlations become unim- 
portant. We experiment with different thin factors (tak- 
ing every step, every lOt/i step and every 50th step and 
we find very little difference in our results. Hence we 
conclude that the sampling of our chains are not greatly 
effected by correlations. 

Lastly, we apply a numerical diagnostic similar to that 
used by Dick et al. [29[ to test the conversion of our 
chains. (This technique is a modification of the Geweke 
diagnostic [30].) We compare the means calculated from 
the first 10% of the chain (after burn-in of 1000) to the 
means calculated from the last 10%. If the chain has con- 
verged to the stationary distribution, then these values 

i i j i j i i T r meani(8i) — mean2(6i) 

should be approximately equal. If ^ ' ±L - LL 

is large, where an is the standard deviation determined 
by the chain for the parameter 6i, then the chain is 
likely to still be drifting. We find that for our chains 

mean l( e,)-~ 2 (e t ) < ^ for Q5% Qf ^ paramcters . The 

remaining parameters are no less than ^P- away from each 
other. Coupling this with the qualitative monitoring of 
the chains described above, we are confident that our 
chains do a good job of reflecting the posterior probabil- 
ity distribution of our model. 



ACKNOWLEDGMENTS 

We thank Matt Auger, Lloyd Knox, and Michael 
Schneider for useful discussions and constructive criti- 
cism. Thanks also to Jason Dick who provided much use- 
ful advice on MCMC. We thank the Tony Tyson group 
for use of their computer cluster, and in particular Perry 
Gee and Hu Zhan for expert advice and computing sup- 
port. Gary Bernstein provided us with Fischer matrices 
suitable for adapting the DETF weak lensing data mod- 
els to our methods, and David Ring and Mark Yashar 
provided additional technical assistance. This work was 
supported by DOE grant DE-FG03-91ER40674 and NSF 
grant AST-0632901. 



[1] D. Huterer and M. S. Turner, Phys. Rev. D60, 081301 

(1999), astro-ph/9808133. 
[2] D. Huterer and M. S. Turner, Phys. Rev. D64, 123527 

(2001), astro-ph/0012510. 
[3] A. Albrecht et al. (2006), astro-ph/0609591. 
[4] A. Albrecht and G. Bernstein, Phys. Rev. D75, 103003 

(2007). 

[5] D. Huterer and H. V. Peiris, Phys. Rev. D75, 083503 

(2007), astro-ph/0610427. 
[6] M. Schneider, L. Knox, H. Zhan, and A. Connolly, As- 

trophys. J. 651, 14 (2006), astro-ph/0606098. 
[7] J. A. Frieman, C. T. Hill, A. Stebbins, and I. Waga, Phys. 

Rev. Lett. 75, 2077 (1995), astro-ph/9505060. 
[8] B. Bozek, A. Abrahamse, A. Albrecht, and M. Barnard 

(2007), arXiv:0712.2884 [astro-ph]. 



[9] M. Barnard, A. Abrahamse, A. Albrecht, B. Bozek, and 

M. Yashar (2007), arXiv:0712.2875 [astro-ph]. 
[10] S. M. Carroll (2001), astro-ph/0107571. 
[11] E. J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. 

Phys. D15, 1753 (2006), hep-th/0603057. 
[12] V. Sahni, Class. Quant. Grav. 19, 3435 (2002), astro- 

ph/0202076. 
[13] R. Bousso (2007), arXiv:0708.4231 [hep-th]. 
[14] P. J. Steinhardt, Phys. Scripta T117, 34 (2005). 
[15] A. Albrecht, AIP Conf. Proc. 957, 3 (2007), 

arXiv:0710.0867 [astro-ph]. 
[16] N. Kaloper and L. Sorbo, JCAP 0604, 007 (2006), astro- 

ph/0511543. 

[17] C. F. Kolda and D. H. Lyth, Phys. Lett. B458, 197 
(1999), hep-ph/9811375. 



11 



[18] S. M. Carroll, Phys. Rev. Lett. 81, 3067 (1998), astro- 
ph/9806099. 

[19] L. J. Hall, Y. Nomura, and S. J. Oliver, Phys. Rev. Lett. 

95, 141302 (2005), astro-ph/0503706. 
[20] M. Dine (2001), hep-th/0107259. 

[21] T. Banks, M. Dine, P. J. Fox, and E. Gorbatov, JCAP 

0306, 001 (2003), hep-th/0303252. 
[22] D. Gamerman, Markov Chain Monte Carlo: Stochastic 

Simulation for Bayesian Inference (London: Chapman 

& Hall, 1997). 

[23] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, 
A. H. Teller, and E. Teller, J. Chem. Phys. 21, 1087 
(1953). 



[24] W. K. Hastings, Biometrika 57, 97 (1970). 

[25] R. R. Caldwell and E. V. Linder, Phys. Rev. Lett. 95, 
141301 (2005), astro-ph/0505494. 

[26] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002). 

[27] N. Christensen, R. Meyer, L. Knox, and B. Luey, Class. 
Quant. Grav. 18, 2677 (2001), astro-ph/0103134. 

[28] B. Carlin and M. K. Cowles, Journal of the American 
Statistical Association 91, 883 (1996). 

[29] J. Dick, L. Knox, and M. Chu, JCAP 0607, 001 (2006), 
astro-ph/0603247. 

[30] E. Geweke, Evaluating the accuracy of sampling based ap- 
proaches to the calculation of posterior moments (1992). 



